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i-Q \ Abstract. The two newly observed hidden-bottom mesons Zj^ (10610) and 

Z^^ (10650) with quantum numbers J^ = 1+ are considered as hadronic molecules 
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Cn . composed of BB* and B*B*, respectively. We give predictions for the widths of 

the strong two-body decays Zj^ — > T{nS) + 7r+ and Z^^'^ — > T{nS) 
Q^ ' phenomenological Lagrangian approach. 
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1. Introduction 



Recently, two hidden-bottom charged meson resonances were observed by the Belle 
Collaboration |T] as two narrow resonance structures in the 7r''=T(n5') {n = 1,2,3) 
and TT'^hfy^mP) {m = 1, 2) mass spectra. They are produced in association with a 
single charged pion in T(5S') decays with the following values of mass and width: 
M[Z+(1061)] = 10608.4 ±2.0 MeV, r[Z+(10610)] = 15.6 ± 2.5 MeV, M[Z^+(10650)] = 
10653.2 ± 1.5 MeV, r[Z^+(10650)] = 14.4 ± 3.2 MeV. Analyses of the charged pion 
angular distributions favor the I^{J^) = l^(l^) quantum numbers of the Z-states [1]. 

Theoretical structure assignments for these hidden-bottom meson resonances were 
proposed immediately after their observation [2]- [10], mainly based on molecular [2] and 
tetra-quark interpretations [9l [10] using the analogy to the charm sector. Also, in [5] 
the new resonances were identified as a hadro-quarkonium system based on the channel 
coupling of light and heavy quarkonia to intermediate open-flavor heavy-light mesons. 

In this paper we analyze the two-body strong decays T(nS')7r+ of Z^^ and Z^~^ 
using a phenomenological Lagrangian approach developed in Refs. [II]-[I6] which is 
based on the compositeness condition [171 [IB]- In particular, in [II]-[l6] recently 
observed unusual hadron states (like D*o{'^317), ^^1(2460), X(3872), y(3940), F(4140), 
Z(4430), Ac(2940), Sc(2800)) were analyzed within the structure assumption as hadronic 
molecules. The compositeness condition implies that the renormalization constant of 
the hadron wave function is set equal to zero or that the hadron exists as a bound 
state of its constituents. It was originally applied to the study of the deuteron as a 
bound state of proton and neutron [17| (see also Ref. [12] for a further application of 
this approach to the case of the deuteron) . Then it was extensively used in low-energy 
hadron phenomenology as the master equation for the treatment of mesons and baryons 
as bound states of light and heavy constituent quarks (see e.g. Refs. [HI [HI [20]). 
By constructing a phenomenological Lagrangian including the couplings of the bound 
state to its constituents and the constituents to other final state particles, we evaluated 
meson-loop diagrams which describe the different decay modes of the molecular states 
(see details in [TT]). 

In the present report we proceed as follows. In Sec. II we briefly review the basic 
ideas of our approach. We now consider the two new resonances Z^ and Z^"^ as the 
two molecular states of BB* and B*B*. Then we proceed to estimate their strong two- 
body decays Z^^ — )■ T{nS) + 7r+ and Z^^ — )■ T{nS) + 7r+ where n = 1,2,3 based on 
an phenomenological Lagrangian approach. In Sec. Ill we present our numerical results 
and a short summary is given in Sec. IV. 

2. Phenomenological Lagrangian approach 

Here we briefly discuss the formalism for the study of the composite (molecular) 
structure of the Z^ and Z^"*" resonances. In the following calculation we adopt the 
spin and parity quantum numbers J^ = 1^ for the two resonances Z^ and Z^,^. We 
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consider these two new charged hidden-bottom meson resonances as a superposition of 
molecular states of BB* and B*B* as 



|Z+(10610)) = -^ 



5*+5° + B*^B^ 

*0\ 



|Z+ (10650)) = \B*+B*''). (1) 

Our approach is based on an interaction Lagrangian describing the coupling of the Z^ 
(or Z^^) to its constituents. The simplest forms of such Lagrangians read 

CzS^) = ^ Mz, Z^{x) I d'y ^zAy') {b{x + y/2)B;{x - y/2) 

+ B;{x + yl2)B{x-y/2)), (2) 

Cz'^x) = ^-^ie,,.pd^z'ax) j dS^z'Sy')B*''{^ + ymB*^{^-y/'^). (3) 

where y is a. relative Jacobi coordinate, g^ and g , are the dimensionless coupling 

constants of Z^ and Z^^ to the molecular BB* and B*B* components, respectively. 
Here ^Zi,{y'^) and $z'(y^) are correlation functions, which describe the distributions 
of the constituent mesons in the bound states. A basic requirement for the choice of 
an explicit form of the correlation function ^niy"^) {H = Z^, Z^^) is that its Fourier 
transform vanishes sufficiently fast in the ultraviolet region of Euclidean space to render 
the Feynman diagrams ultraviolet finite. We adopt a Gaussian form for the correlation 
function. The Fourier transform of this vertex function is given by 

^Hipl/A') = expi-pl/A') , (4) 

where pe is the Euclidean Jacobi momentum. A is a size parameter characterizing the 
distribution of the two constituent mesons in the Z^ and Z^^'^ systems, which also leads 
to a regularization of the ultraviolet divergences in the Feynman diagrams. From our 
previous analyses of the strong two-body decays of the X, Y, Z meson resonances and of 
the Ac(2940) and Sc(2880) baryon states we deduced a value of A ~ 1 GeV [15j. For a 
very loosely bound system like the X(3872) a size parameter of A ~ 0.5 GeV |!l4j is more 
suitable. The coupling constants g^ and g , are then determined by the compositeness 

*> b 

condition [171 dSl [T^l [13 [TT] . It implies that the renormalization constant of the hadron 
wave function is set equal to zero with: 

Zh = 1-^h{MJj) = 0. (5) 

Here, T^'^j is the derivative of the transverse part of the mass operator E^^ of the 
molecular states (see Fig.l), which is defined as 

S^^(p) = ^rs^(p) + ^Si(p), gT = 9"'-^. (6) 

The compositeness condition (|5]) gives a constraint on the choice of the free parameter A. 
Analytical expressions for the couplings g and g , are given in Appendix A. In the 

b Zjj 

calculation the masses of Zb and Z^ are expressed in terms of the constituent masses 
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Figure 1. Mass operators of Z^ and Z^^. 
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Figure 2. Two-body decays Z+ -^ T{nS) + tt and Z(,+ -> T{nS) + n 



(7) 



and the binding energy e (a variable quantity in our calculations): 

Mz, = MB + MB*-e, Mz^ = 2Mb* - e . 

Here we assume bound states for the Zb and Z^. 

In the calculation of the two-body decays Z^^ — )• T(r7,5')-|-7r+ and Zj^^ — t- T{nS) + n^ 
we include the direct four-point interactions for the BB*Ttt'^ and B*B*T7i^ vertices. 
The respective phenomenological Lagrangians take the form 

CBB*Tn{x) =g^^,^T^{x)B*^{x)7i{x)-rB{x) + Ii.c., (8) 

CB*B*rM = te,,^^(^g^^^,^T''{x)B*%x)d''n{x) ■ tB*''{x) 

+ f^a'^T^{x)B*^{x)n{x)-TB*^ 



X, 



(9) 



The four-particle coupling constants defined in Eqs. ([8]) and iQ are effective, which 
also include off-shell effects. Such couplings are obviously different from the strong 
couplings of molecular states to their constituents which model their composite structure 
via vertex function distibutions. We use effective Lagrangians (using both SU(4) and 
SU(5) classification schemes) developed by Ko and collaborators [21] which worked 
phenomenologically quite successfully. Corresponding diagrams contributing to the 
Zf^ — )■ T{nS) + TT"*" and Z^^ — )■ T{nS) + n^ processes are shown in Fig.2. 

Among the three couplings (7^^,^^, 9^*^,^^ and fg*g*^^ we have the relations [2T] 

_ 9tBB 9b*Btt „ _ f _ dsB'TTT 



The hadronic couplings g^gg and g 



^B*B-it{x) 
^TBBix) 



TBB 

9f 



are defined as 



2 a/Ms Ms* 



(10) 



B*Jx)id>'n{x)rB{x)+li.c. 



V2 
g^gg T^{x)B{xWB{x) + H.c. 



:ii) 
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The coupling constant g^,g.gg is given by 

/T(n5) 

where fr(nS) is determined from the leptonic decays of the T{nS) states as 

rfT(nS)^e+e-)=^^#^, (13) 

\ ^ ' J 27 Mrins) ^ ' 

where a^^ = 1/137.036 is the fine-structure constant. The relation (fT2|) is analogue of 
the p-meson universality 

«., - ^ - f (14) 

extended to the bottom sector in Ref. [2Tj, where Qp^ is the p — )■ 7 transition 
coupling. Here For the last couplings we get /t(i5) = 715.2 MeV, /t{25) = 497.5 MeV, 
/^(3g) = 430.2 MeV, where we used the mass values Mx(is,2s,3s) = 9460.30 ± 0.26 MeV, 
10023.26 ±0.31 MeV and 10355.2 ±0.5 MeV as well as the results for the leptonic decay 
widths of the T{nS) states 

r h^lS) -^ e+e-\ = 1.340 ± 0.018 keV , 
r (t{2S) -^ e+e") = 0.612 ± 0.011 keV , 

r (t{3S) ^ e+e~) = 0.443 ± 0.008 keV . (15) 

Note that we explicitly take into account the My(^nS) dependence of the fr{nS) and g^gg 
couplings. 

The coupling g^g*^ can be related to the effective coupling constant g = 0.44 ± 

0.01 
0.00 



0.03^'^ogo determined in a lattice calculation [22]. The relation is 

gB*Bn = J- VMbMb* ~ 35.34, (16) 

where fj^ ^ 132 MeV is the pion decay constant. We should stress that the 
phenomenological Lagrangians developed in Refs. [21] were successfully applied to 
different aspects of heavy flavor physics. 

3. Numerical results 

With the phenomenological Lagrangians introduced and discussed in Sec. II one can 
proceed to determine the widths of the two-body decays Z^ — )■ T{ns) + vr"*" and 
Zfj'^ — )■ T(ras) ± TT"*" with n = 1,2,3 (see corresponding diagrams in Fig.2). The 
corresponding decay widths are given by: 

p 9ZtT(nS)TT ^1/2^ „ ^2 »/f2 iy^2\ 



(17) 



rz-.x(..v.-fe^AV^(M|,M^(„,,M^) 



,TinS)n+ IQt,Mz 
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where X{x, y, z) = x^+y'^+z^—2xy—2xz—2yz is the Kallen function. The decay couphng 
constants gz^T{nS)-n and gz'^T{nS)n involve the products gz^T{nS)7r = 9z, S'ss-rcns). ^^ 
and gz'jinS)7r = 9^, 9b*b*i:(^s)^ ^z^ J2 where the loop integrals Ji and J2 are given in 
Appendix A. 

For our numerical evaluation hadron masses are taken from the compilation of the 
PDG [23]. The only free parameter of our calculation is the dimensional parameter A 
entering in the correlation function of Eq. (jl]). As mentioned before, the parameter 
A describes the distributions of BB* and B*B* in the Z^ and Z^^ bound state 
systems, respectively. Tables I and II contain our estimates for the decay widths of 
Zf^ —7- T{nS) + 7r+ and Zj^^ — t- T{nS) + 7r+. We also indicate the values for the 
couplings g^ and g , as determined from the compositeness condition. We find that 
the data on the strong T{nS)7!' decays of Z^ and Z^ can be qualitatively described 
for A ~ 0.5 GeV. This value is close to the one used for the molecular system of the 
X(3872)[T^. This also means that the Zf, and Z^ states are considered as extended 
molecular states. For transparency we present the results for several choices of the size 
parameter A = 0.4, 0.45, 0.5, 0.55 GeV and the binding energy e. Note that an increase 
of A leads to a larger decay width. Although the dependence of the decay widths on the 
binding energy is rather moderate, a quantitative prediction for these decays strongly 
depends on the size of the system. The decay pattern of Z^ and Z^ to T{nS) + 7r+, 
that is the relative importance of the T{nS) + 7r+ decay channels for n = 1,2,3, can 
be reproduced and could give further support for the molecular interpretation of these 
states. One can see from Tables I and II that the rates are increased by a factor 2-2.5 
which means that the amplitudes are roughly increased by a factor 1.4-1.6. Latter value 
is consistent with a growth of the cutoff parameter of 0.55/0.4 = 1.375. In this respect 
the calculation is consistent and we do not see any disagreement. It is also clear that 
any model based on cutoff regularization has cutoff-dependent result. Here the cutoff is 
related to the size of hadronic molecular compound. As done in our previous analyses 
of other heavy hadron molecules data help to do a fine tuning of the cutoff parameter 
which is specific for a particular molecular state. There is no a universal cutoff for the 
whole tower of possible hadronic molecular states. Our previous analyses of molecular 
states consisting of two heavy mesons (like the X(3872) state) indicate that the cutoff 
parameter in the vertex functions of Z^ and Z^ states should be around 0.5 GeV. This 
is the reason why we do predictions for the valus of A close to 0.5 GeV and are also 
waiting for more precise data with smaller error bars. Please also note that the ratio of 
rates is less dependent on the cutoff and hence a stabile prediction of the model. Since 
data are given for the absolute rates we also choose to use this presentation for our 
results and not the relative rates. 
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Table I. Z^ — )> T(nS') + 7r+ decay properties. 



e(MeV) 


9z^, 


ri5(MeV) 


r25(MeV) 


r35(MeV) 


1 


3.3, 3.4, 3.6, 3.7 


11.6, 16.4, 22.3, 29.5 


13.7, 19.4, 26.4, 34.9 


7.2, 10.2, 13.9, 18.4 


5 


4.0, 4.0, 4.1, 4.2 


11.0, 16.5, 22.4, 29.5 


13.0, 19.4, 26.3, 34.7 


6.7, 10.1, 13.6, 18.0 


10 


5.0, 4.9, 4.8, 4.8 


11.7, 16.9, 22.7, 29.7 


13.7, 19.8, 26.6, 34.8 


7.0, 10.1, 13.5, 17.6 


20 


7.2, 6.6, 6.3, 6.0 


13.3, 18.3, 24.0, 30.8 


15.4, 21.2, 27.8, 35.7 


7.4, 10.2, 13.4, 17.3 


30 


9.4, 8.5, 7.9, 7.4 


14.5, 19.6, 25.4, 32.3 


16.7, 22.5, 29.2, 37.1 


7.7, 10.3, 13.4, 17.0 


40 


11.7, 10.4, 9.5, 8.8 


15.4, 20.4, 26.7, 33.7 


17.5, 23.2, 30.3, 38.3 


7.5, 10.0, 13.1, 16.5 


50 


13.9, 12.3, 11.1, 10.2 


15.9, 21.2, 27.6, 34.9 


17.8, 23.9, 31.1, 39.2 


7.1, 9.6, 12.5, 15.7 


Exp. 




22.9 ±7.3 ±2 


21.1 ±41^ 


12.2 ±1.7 ±4 


Table II. Zj^ — )> T{nS) + 7r+ decay properties. 


e(MeV) 


gz'. 


ri5(MeV) 


r25(MeV) 


r35(MeV) 


1 


3.2, 3.4, 3.6, 3.7 


12.0, 16.9, 23.0, 30.4 


14.7, 20.8, 28.3, 37.4 


9.0, 12.8, 17.4, 23.0 


5 


3.9, 4.0, 4.1, 4.2 


12.1, 17.0, 23.0, 30.3 


14.9, 20.8, 28.2, 37.2 


9.0, 12.6, 17.1, 22.6 


10 


5.0, 4.8, 4.8, 4.7 


12.7, 17.4, 23.4, 30.6 


15.4, 21.3, 28.5, 37.3 


9.2, 12.7, 17.1, 22.3 


20 


7.2, 6.6, 6.3, 6.0 


14.0, 18.8, 24.6, 31.7 


16.9, 22.7, 29.8, 39.3 


9.8, 13.2, 17.3, 22.3 


30 


9.4, 8.5, 7.8, 7.3 


15.1, 20.1, 26.1, 33.1 


18.1, 24.1, 31.3, 39.7 


10.2, 13.5, 17.6, 22.3 


40 


11.7, 10.4, 9.4, 8.7 


15.8, 21.1, 27.3, 34.5 


18.9, 25.1, 32.4, 41.0 


10.2, 13.6, 17.6, 22.2 


50 


13.9, 12.3, 11.0, 10.1 


16.2, 21.7, 28.2, 35.6 


19.1, 25.5, 33.2, 41.9 


9.9, 13.3, 17.3, 21.8 


Exp. 




12 ±10 ±3 


16.4 ±3.6^^ 


10.9 ±2.6+^ 



4. Summary 

To summarize, we have pursued a hadronic molecular interpretation for the two recently 
observed hidden-bottom charged mesons Z^ and Z^^ . In our calculation we have used 
the spin-parity assignment J^ = 1+ for the two resonances, which is currently favored 
by the experimental decay distributions. We have studied the consequences for their 
two-body decays Z^ — )■ T{nS) +7r+ and Z^^ — )■ T{nS) + 7r"^ within a phenomeno logical 
Lagrangian approach. The calculated results for the decay widths (see Tables I and 
11) are of the order of MeV and for the most part qualitatively consistent with the 
numbers deduced by the Belle collaboration. Especially the experimental results for 
Z^(Zj,^) — 7> T(nS') + 7r+ can be reproduced taking a value for the free size parameter A 
near 0.5 GeV. 
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Appendix A. Coupling constants and structure integrals 

The expressions for the coupling constants g^ ,g , and structure integrals Ji, J2 are 

b Z^ 

r, CO 

_2 M^^ f daida2 , , „ \ A , ^ 



"b 



9~ = TTTT^TT^ / A^ ("12 + 2aia2) 1 + 



327r2A2 j A? ^ '' ' ' " V 2M|.Ai 







r M|,ai + M|a2 Mf ai2 + 2aia2l ,.,, 

"^^n A^ +2A^ A, I' (^-'^ 

T\./r2 °° / \ 

_2 -^b /" (iai(ia2 / A^ ai2 + 2aia2 \ A , A^ 

9 71 ~ 1C_2 A2 / a1 772 '" 



167r2A2 y A? [Ml, 2Ai / V M|.Ai 







1 /" daida2 f A^ 

-'l = 777- / ^7 — 1 



57r2 y A2 V 2MI.A2 



X exp I ^^ + jj^ X (^ ' (A-3) 

00 
1 /" (ia;i(ia;2 / A^ 

<^2 = :: — 7: / — rT5 — 1 + 



7r2 y Ai V MI.A2 



r2 



M|.ai2 ^z^ ai2 + 4aia2 

^V^ 4A^ A^ 



xexp^-i:l^f^ + ^ -7— M , (A.4) 



where 



Ai = 2 + ai2 , A2 = 1 + q;i2 , a;i2 = cti + 02 • (A. 5) 
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